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ABSTRACT: 

The requirements of a track -while-scan radar data 
processing scheme are stated from the point of view of the 
evader in a pursuit-evasion game. The rules of the game, 
determined in part by the nature of the pursuer's weapons, 
are such that the evader must be able to discriminate re- 
liably between straight-line and maneuvering pursuer 
motion. A suggested method for such discrimination is 
tested by simulation. The method employs bias-sensitive 
maneuver detection and gain-adaptive discrete Kalman 
filtering. Also tested is a smoothing scheme for esta- 
blishing long-term trends in a pursuer's maneuvering track. 
The outcome of both tests indicate that the suggested pro- 
cessing methods may be useful in the formulation of fire 
control policies for destruction of maneuvering targets. 
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Introduction 
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Consider a pursuit-evasion problem where the pursuer is 
a high-speed, highly maneuverable object and the evader is 
slower and less maneuverable. Both pursuer and evader possess 
weapons capable of destroying their respective adversary. For 
purposes of illustration, we shall assign the role of pursuer 
to a PT-type boat, and that of evader to a ship the size of a 
destroyer or larger. Both are presumed to possess radar for 
purposes of relative position location. 

This paper concentrates upon the tracking methods employed 
by the evader, whose tactical decisions will be conditioned in 
larger measure upon hopefully reliable estimates of the pursuer's 
current and near-future intent. A significant measure of this 
intent is the nature of the pursuer's present motion, i. e., is 

■v 

he maneuvering or is he proceeding in a straight line? The 
significance of maneuver information lies in the fact that the 
pursuer, for purposes of proper gyro alignment, must proceed 
in a straight line for a given time interval just prior to 
launching a torpedo. Given a tracking scheme with reliable 
maneuver detection, the evader may thwart the pursuer's attempted 
launch by firing projectiles in such a manner as to destroy the 
pursuer or at least force him to maneuver and thus abort the 
launch. 



The evader's track-while-scan/maneuver detection requirements 



We shall assign two conditions of motion to the pursuer; 
he is either maneuvering or moving in a straight line. The 
maneuver detector will at any given time, then, fall into one 
of the following four possibilities: 



Detector says : 
A straight 
B maneuvering 
C straight 
D maneuvering 



Pursuer is : 

straight 

maneuvering 

maneuvering 

straight 
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Of these possibilities, A and B are obviously to be desired. 

Of the two error conditions C and D, D is by far the more 
dangerous because it allows the pursuer to carry out a weapon 
launch without the harassment of the specific anti-launch 
tactics that would otherwise be carried out by the evader. 

Error condition C, on the other hand, is an error in the con- 
servative direction, and may result at worst in wasted ammu- 
nition. Minimization of D-type errors is therefore stated as 
a requirement. 

For purposes of further discussion, two descriptors are 
now defined. The descriptor "short-term” shall apply to esti- 
mates and predictions of pursuer position and velocity no more 
than one sample interval ahead of the current time. The des- 
scriptor "long-term" shall apply to estimates and predictions 
of the same quantities many samples ahead of the current time, 
where "many" is understood to be related to the flight time of 
ballistic projectiles fired by the evader. 

Common sense tells us that it would be foolish to attempt 
a classical fire-control solution (i. e., one in which the 
target is assumed to proceed from a known position with constant 
speed and heading) when the pursuer is maneuvering. This is not 
to say that the evader should, under these conditions, refrain 
from firing shells; it does say that he needs long-term pre- 
diction capability of a markedly different nature than that 
which . would be obtained by simple extension of short-term 
predictions. Long-term prediction under maneuvering conditions 
will require considerable smoothing of data from the immediate 
and extended past, to the end that gross trends in the pursuer's 
track are established. Once obtained, such trends can form the 
basis for patterns of fire covering areas of ocean in which it 
is most probable that the pursuer will be after shell flight 
time has elapsed. 

The requirement for short-term tracking is two-fold. The 
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maneuver sensor will base its decisions on the difference 
between the pursuer's radar returns for a series of observa- 
tions and the corresponding short-term predicted positions for 
that series. Secondly, the short-term tracker must be such 
that after tracking a non-maneuvering target for a given length 
of time, its velocity estimates must converge to such an accuracy 
as to allow extension to long-term position prediction for the 
gun fire control solution. 

Fig. 1 summarizes, in block diagram form, the various 
functions to be accomplished in the track -while-scan processing 
described above. The following sections will describe the 
details of a particular processing scheme that was simulated 
and tested on the digital computer. 

A suggested track-while-scan processor 

The gross features of the suggested processor are given 
as follows: 

Short-term tracker -- a Kalman filter (sequential 

estimator) with provision for backsliding in the gain 
schedule consistent with maneuver information. 

Long-term tracker — a least-squares fit of a straight 
line, constant velocity track to at least 40 seconds 
worth of immediate past data 

Maneuver sensor -- an indication is given that a maneuver 
may be in progress if 

a) two consecutive differences between predicted 
and radar x or y locations exceed a gate of 
given size 

and b) both differences are of the same sign. 

(Detailed descriptions of the processing are included in the 
following section on simulation) . 
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EVADER GUN-FIRE COMMANDS 



Fig* 1 Block diagram of evader's track— while— scan processing 
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A simulation for evaluation of track-while-scan processing 



A digital computer simulation program was written to in- 
vestigate and evaluate the processing methods outlined above. 

Pursuer Tracks 

Two standard pursuer tracks were used in the evaluation. 
Both tracks were initiated at a range of 15,000 yards and an 
azimuth of 45°. The pursuer maintains a speed of 30 knots in 
both tracks. In track A, the pursuer proceeds toward the eva- 
der (at the origin and stationary) in straight-line motion for 
the entire 80 second track history. In track B, the pursuer 
starts out as in Track A, but commences a zig-zag maneuver at 
20 seconds. The zig-zag continues for 40 seconds, whereupon 
the track resumes straight line motion to 80 seconds total 
time. A and B are shown, with broken axes, in Figure 2. 

Evader's Radar 

The evader's radar is assumed to have a one-second 
data rate. Radar measurement errors are assumed to be 
gaussion, zero-mean, white sequences with standard deviations 
of 4 yards in range and 2 milliradians in azimuth. The evader's 
digital processor is assumed to accept noisy range information 
with 10-yard quantization levels, and noisy azimuth information 
with 13-bit encoder accuracy. These quantizations are included 
in the simulation. 



Short-term Tracker 

The short-term tracker is basically a standard Kalman 
filter based on a second-order, pure inertia target dynamics 
description in uncoupled x,y coordinates. In pre-computing the 
filter gain schedule, the cartesian measurement error covariance 
matrix was made diagonal and constant indicating that we assume 
the (x,y) measurement errors to be additive and independent. 
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Fig. 2 Standard straight-line and aig-aag pursuer courses 
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This is, of course, not the case, since coordinate transforma- 
tion from polar (r,0) to cartesian (x,y) involves the polar 
errors multiplicitavely and nonlinearly in the cartesian errors. 
Furthermore, cartesian error components are space dependent, so 
that ctx and ay are non-constant. To fully account for all these 
effects of the coordinate transformation upon cartesian measure- 
ment error components would require on-line computation of the 
error covariance and of filter gain values, at considerable cost 
in computation time. Filter performance with a gain schedule 
based on constant and independent cartesian error statistics 
has proved quite satisfactory; the results are felt to justify 
the assumptions. At worst, if a given tracking situation should 
result in a large spread of pursuer-evader relative range, two 
or more complete gain schedules could be stored and used for the 
appropriate range levels. 

The filter prediction-correction equations are identical 
for both x and y, and are as follows: 

Prediction jy k _j = k . 1 + -T (1) 

*kl k-1 = *k-ll k-1 (2) 

Correction Sy k - ^ k . x + 9** • (* k - ^ k _0 (3) 

*k| k = *k| k-1 + g2 k C z k ~ *k| k-l) (4) 

where z. is the radar observation in x or y. 

K 

The gain schedule (gl,g2)^and the equations and constants 
used in its calculation are given in Appendix A. It should be 
noted here that a null random excitation covariance was used in 
the calculation of the gain schedule. This places the necessity 
of following pursuer acceleration with the maneuver detector. 
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and allows the short-term filter to do a more accurate job of 
estimation when the pursuer is moving in straight-line motion. 

One consequence of using a null excitation covariance is the 
convergence of filter gains to zero value. To allow the actual 
schedule to so converge, without truncation, would be unwise 
because long-term drift tendencies in a nominally straight 
pursuer track could lead to highly undesirable diverging biases 
in the estimation. The use of a drift detector has been suggest- 
ed and is certainly feasible. The present study, however, does 
not include drift detection. 

Built into the short-term tracker is a backsliding feature 
which causes a re-processing of the most recent three pieces 
of radar information. The rationale for backslide processing 
is this; the maneuver detector indicates positively when it 
senses a sizable estimation bias building up over two consecu- 
tive observations. The ultimate effect of that bias, and the 
time required to recover from it once it is recognized, can be 
reduced by going back to where the bias most likely started 
accumulating, namely two sample intervals back. 

If the radar data collected two samples back is then given 
more credence in the re-processing, the effect will be to reduce 
the bias. This increase in the weighting of back radar data 
can be accomplished by going back to a very early point in the 
gain schedule, where respective gains are associated with rela- 
tively infirm tracks. In maneuvering situation, the track 
becomes degraded. 

Let us be explicit about the nature of the re-processing. 

til 

Let us say that we have just received the k radar observa- 
tion, we have compared it with the conditional estimate or 
prediction of what the observation should have been, and we 
have decided on the basis of this and the previous comparison 
that a maneuver is in progress. The following equations are 
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solved, in the order given. 



*k-2| k-2 *k-2| k-3 + gl n * ( xR k-2 *k-2| k- 3 ) 

K-2\ k-2 = \-2\ k-3 + g2 n * ( xR k-2 “ *k- 2 | k- 3 ) 
^k-ll k-2 = *k-2| k-2 + ^k-21 k-2 • T 



^k-ll k- 


-2 = K-2\ 


k-2 
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K\ k-i + 


g2 n+2 


( xR k - 


*kl k-i) 







This set of equations differs from the basic prediction 
correction set (1) through (4) in the gain index. In equations 
(5), we note that the re-processing starts with a gain schedule 
index of n, where n ^ k-2. It was found by trial and error in 
the simulation that n = 1 was the most satisfactory starting 
point in the gain schedule for re-processing. After the execu- 
tion of equations (5), routine processing continues, with the 
gain schedule index proceeding to n+3, n+4, etc., until a 
maneuver indication is once again received. 
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The Maneuver Detector 



On the basis of a 100-member ensemble, raw radar carte- 
sian errors for the standard tracks of Figure 2 were found to 
be 20 yards standard deviation in both x and y. By trial in 
the simulation, a maneuver detector gate of ± 40 yards was 
found to be most satisfactory. This is equivalent to ± 2cr on 
the raw radar. The detector was made to give a re-process/ 
backslide command when, for two consecutive observations, the 
predicted and observed cartesian positions differed by more 
than 40 yards, with the differences of the same sign. In 
equation form. 



and both differences have the same sign, then a maneuver indi- 
cation is given. 



This test consisted of processing an ensemble of 100 runs 
of standard pursuer track A of Fig. 2, i. e., straight-line, 
constant velocity motion. Each run differed from the others 
in the measurement noise sequence used in the generation of 
simulated radar observations. The test was conducted to deter- 
mine two things; 1) how often would the maneuver detector give 
a positive (in this case, false) indication, and 2) what sort 
of runs errors in position and velocity could be expected from 
the short-term tracker operating on an extended straight-line 
tracks. The first result was obtained simply by establishing 
a counter to count the number of re-processing commands gene- 
rated by the detector. This result was 53 commands. Expressed 
as a percentage of the total number of observations processed. 



if 




( 6 ) 



and 




Non-Maneuvering Test 
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which was 100x80 = 8000, the detector gave false indications 

on 0.665% of the observations. This percentage could be 
reduced even further by requiring two or more consecutive 
re-process commands to occur before the track is interpreted 
as maneuvering. This does not suggest that the indicated 
re-processing be suppressed, but that several re-process 
commands in sequence be required before stating that the 
pursuer is maneuvering. 

Figures 3 and 4 are plots of the ensemble rms error in 
x and x, respectively, for estimates produced by the short- 
term tracker. Embedded in these plots are the effects of the 
0.665% false maneuver indications as discussed above, i. e., 
re-processing was done where indicated. The steady state 
values of these rms errors are approximately 6 yards and 1 
yard per second, respectively. Similar plots were obtained 
for the y coordinate. 

Maneuvering Test 

This test consisted of processing an ensemble of 100 runs 
of standard pursuer track B of Figure 2, i. e., a combination 
of straight and zig-zag motion at constant speed. As in the 
non-maneuvering test, each run differed from all others in the 
measurement noise sequence used in the generation of simulated 
radar observations. The test was conducted to determine the 
effectiveness of the short-term tracker/maneuver detector 
against a maneuvering pursuer, and to observe the trend-esti- 
mating capability of the long-term, least-squares tracker. 

Figure 5 is a plot of the standard maneuvering track 
upon which are superimposed vertical bars whose lengths are 
proportioned to the number of times, in the ensemble of 100, 
that re-processing was dictated by the maneuver sensor at that 
particular point in the true track. It may be noted that the 
maneuver sensor quite rapidly re-establishes a non-maneuver 
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RMS error in estimation of X for an ensemble of tracks (A) 
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RMS error in estimation of X for an ensemble of tracks A 



Scale: one re-process command 
millimeter of vertical 
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Fig. 5 Re-process command history for track B, 
ensemble of 100 trials 
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indication after 60 seconds, during the tail-end straight por- 
tion of track B. 

Figure 6 and 7 are plots of the ensemble rms error in x 
and x, respectively, for estimates produced by the short-term 
tracker. These plots display quite vividly the effect of the 
maneuver on the short-term tracking; for the first twenty 
seconds, they are identical to Figures 3 and 4. The zig-zag 
portion of the track extends over the time interval from twenty 
to sixty seconds, during which time the rms error in position 
rises to a level somewhat below that of the raw radar data. 
After the track resumes straight-line motion at sixty seconds, 
the error levels in both position and velocity resume their 
downward trend toward the values achieved in the steady-state 
of Figures 3 and 4. Had the short-term tracker not been 
adaptive to maneuvers in the sense of gain schedule backslid- 
ing, the true and filtered tracks would have diverged. 



It is obvious from Figure 7 that velocity estimation suffers 
considerable inaccuracy during the maneuver portion of the track. 
This degradation is acceptable, since these maneuvering velocity 
estimates are not to be used in establishing long-term position 
predictions for fire control purposes. Such fire would be 
futile in the light of the evader's knowledge that the pursuer 
is maneuvering. Only when the pursuer is known not to be 
maneuvering are the short-term tracker's velocity estimates 
used for fire control position prediction. 

The Long-Term Tracker 

During those periods in which the pursuer is maneuvering, 
the evader may wish to lay down a pattern of fire based not 
upon point-by-point predictions of pursuer position, but upon 
areas or zones through which it is highly likely that the 
pursuer, over a given time interval, will pass. The location 
of these zones may in some measure be aided by determining 
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BMS error in estimation of X for an ensemble of tracks B 



RMS error, yards/second 
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long-term trends in the pursuer's motion. The feasibility of 
processing radar data so as to provide this trend information 
was tested by applying least squares curve fitting techniques 
to the same ensemble of maneuvering tracks described above. 

The processing was done as follows. As soon as the first fourty 
radar observations were taken and stored, a straight line, con- 
stant velocity, least squares track was fitted to the observa- 
tions. Thereafter, as a new radar observation became available, 
the most remote data were dropped, thereby keeping the observa- 
tion list of size fourty and effecting a sliding fit over the 
last fourty pieces of data. 

The processing algorithm for such a fit is very straight- 
forward. Straight line pursuer motion may be described (for 
the x cartesian coordinate) as 

x ( t ) = x (0 ) + v x • t (7) 

If we devote the raw cartesian observations at t = t. as x., 

1 l 

the difference between the observation and the fitted line 
will be 

c. = x. - x(t.) = x. - x (0) - v • t. (8) 

1111X1 v ' 

If we sum the squares of all such differences over our list of 
fourty observations, we have 




To minimize this sum of squared differences, 

3 e 5c 

■§^ToT = '§v - = 0 < 10 ) 

' X 
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Satisfaction of equation (10) yields 



40 40 40 

r ,x(0) + 5 ,v • t. = 5 ,x. 
1=1 ' 1=1 X 1 1=1 1 



40 40 2 40 

V n x(0) • t. + £ ,v ♦ t. = £ ,x. • t. 

i=l ' ' 1 1=1 x 1 1=1 1 1 



( 11 ) 



For purposes of convenience, the time reference is selected so 
that t^ = - 39 seconds and t^Q = 0 seconds, with the result 

that x(0) is the least squares estimate of position "now", i.e.,at the 
most recent end of the forty-observation list. The resulting 
solution for present position and velocity is 



x(0) 



40 40 

20,540 ? . x. + 780 £ _ x. • t. 
i=l i i=l i i 

213,200 



40 40 

40 ? , x. • t. + 780 ? . x. 

1=1 i i 1=1 l 

v = 

x 213,200 



(12) 



The least squares positions and velocities obtained through 
equations twelve were averaged over the ensemble for times 41 
through 80 seconds. The x and y velocities are plotted in Figures 
8 and 9, respectively. 

Figure 10 is another copy of the standard maneuvering track 
( B) upon which are superimposed several representative position 
and velocity estimates from the long-term tracker. The origin 
of each arrow is placed at the least-squares position for that 
point in time. The line from that point to the true track 
indicates where the pursuer actually is at this time. The head 
of the arrow rests upon that point at which the track would be 
in twenty seconds, based upon least-squares velocity for position 
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Fig. 9 Mean I from sliding least-squares 
fit, for ensemble of tracks B 
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prediction. The results shown in Figure 10 may be interpreted 
as follows. The first output of the long-term tracker, at 
time 41 seconds, shows the considerable influence of the ini- 
tial straight-line portion of the track. As time progresses 
beyond 41 seconds and much of the straight-line data is drop- 
ped off the stale end of the observation list, the long-term 
predictions begin to align most favorably with the directional 
trend of the track's zig-zag section. This alignment is 
especially evident in the time interval from 55 to 65 seconds. 
For times beyond 65 seconds, the long-term tracker is producing 
obviously poor predictions, but these are of little importance 
because by this time, the short-term tracker and maneuver 
detector have indicated that the pursuer is once again moving 
with straight-line motion, with full tracking responsibility 
shifted to the short-term tracker. 

It should be noted that the structure of the long-term 
tracker used in this study is by no means suggested as the 
"only way to go". There is nothing sacred about a straight- 
line fit; higher-order polynomial curve fitting may provide 
even better results. Nor is fourty seconds duration for the 
observation list necessarily the best figure. The point we 
hope to have illustrated is that in the context of this parti- 
cular set of tactical circumstances, for reasons discussed 
above, it may be highly advantageous to run two parallel radar 
data processors. 

Conclusions 

In a pursuit-evasion game where pursuer straight-line 
motion constitutes the greatest threat to the evader, it is 
essential that the evader have the capability of discriminating 
between straight-line and maneuvering pursuer motion. The 
combination of a maneuver detector and a Kalman filter whose 
gain schedule is made adaptive with respect to detected maneu- 
vers has been shown to be effective both for maneuver tracking 
and for meeting relatively high accuracy specifications for 
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steady-state tracking of straight-line motion. 

It has been shown, also, that parallel processing of 
the evader's radar returns by least squares curve fitting 
may be quite helpful in establishing long-term trends of 
the maneuvering phases of the track. These trends might 
conceivably be used in forming a pattern of fire against 
a maneuvering target. 
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Appendix A 



The recursive Kalman filter equations from which the gain sche- 
dule was computed are given as follows. The equations apply 
to x and y axis estimation alike. 



gl k 


T 

= PH 


HP H T + R 


_ g2 k_ 


k| k-1 


k| k-1 



P k|k = (I - G k H) P k|k-1 



? k+l| k 



= S> p $ + Q 



where the transition matrix $ is given by 



$ = 



1 

0 



T 

1 



the measurement matrix H by 



H = 



t <0 



and the excitation covariance matrix Q by 



(A-2) 



(A-3) 



a - [o] 

The measurement error covariance was chosen as 



(A-4) 



R = 




(A-5) 



on the basis of experimental mean square values of quantized 
radar position minus true position. 



The recursive computations were started with the follow- 
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ing initial condition. 



P l| 0 



25(10 6 ) 0 

0 300 



(A-6) 



The initial filter state vector for each member of the ensemble 
was given the value 

~ x 

x 

The first fifteen seconds of the gain schedule, as computed from 
equations (A-l) through (A-6), is shown in Table A-l. 



ll 0 



0-0 

0-0 



(A-7) 
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k 


gl 


g2 


1 


0.99 


0.00 


2 


0.64 


0.27 


3 


0.63 


0.30 


4 


0.60 


0.24 


5 


0.55 


0.18 


6 


0.49 


0.13 


7 


0.45 


0.10 


8 


0.41 


0.08 


9 


0.37 


0.07 


10 


0.34 


0.05 


11 


0.31 


0.05 


12 


0.29 


0.04 


13 


0.27 


0.03 


14 


0.25 


0.03 


15 


0.24 


0.02 


Table 4-1 


The first fifteen se< 



the short-term tracker gain schedule 
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Appendix B 



Figure B-l is an overall flow chart for the tracking simula- 
tion program employed in the study. 
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